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INTRODUCTION 


A version of the Woodward Carmichael panel pressure computer analysis tech 
niquc used at the National Aeronautics and Space Administration Drydcn Flight 
Research Center (NASA DFRC) has been applied to several diverse aerodynamic 
studies. The purpose of this paper is to document work performed with this com 
puter program. 

The program calculates the potential, inviscid, steady state flow solution for a 
given wing and body combination in subsonic or supersonic flow. The program was 
developed at the NASA Ames Research Center using the Woodward technique. A 
description of this computational method is included in references 1 and 2. The 
latest version of the program has been in use id DFRC since 1975 and is commonly 
referred to as Wing Hotly. 

Experience gained in the use of this program includes Indh its intended appli 
cations as well as some unusual applications This report presents the methods and 
success of using Wing Body to determine aircraft stability and control characteris 
tics, potential (low field characteristics . air launch dynamics, and wake vortex up 
set loads. The limitations and range of application of tin* program are also discussed. 

The greatest advantage of the Wing body version over most existing aerodynamic 
analytical techniques is that it is user oriented. Preparation of computer input is 
simplified by the use of manv default values and an automatic paneling feature. The 
input scheme is also very versatile, allowing nonconvent ional uses of the program. 

The Wing body analytical techniques for determining aircraft stability and con 
trol derivatives are shown to be valid. The ef'eets of design changes can be closely 
predicted prior to wind tunnel and flight testing. Other analytical procedures, such 
as the determination of launch dynamics or wake vortex induced loads, provide 
results which cannot be easily obtained by experiment. 


I'lio analy ats presented m lilts paper was performed l*y mi engineering student 
trainee at PFRi* ns purl of llio cooper at tv o tr Minot' program . 

SYMMOI S 

riio body axis s\ stem assumed throughout this report is consistent with that 
used by Wing body. \ is positive in llio nfl tlirivlion, V is positivo to the right, 
Mini / is positivo up 

A pMMOl area, ITl* 

b roforonoo spun . m 

r, lift coefficient 

li 

rolling moment coefficient 

pitching momont coefficient 

t v , norms l force coefficient 

C' n yawing momont oKflticlont 

l’ proMMuro coefficient 

P 

t’y Stilt’ ft'l't't* Coefficient 

e roforonoo chord , m 

e^ v , chord of M 59 nit* plant*, span station -OS. m 

K number of panels 

M Maoli nuinbor 

n variable 

S roforonoo ptnnform area, m* 

V airspeed, m s*v 

V vortical velocity component . in see 


J 



\i' loi>ditudinal dtstanoo from panol oontronl to oontor ot d«*«' u> . m 

NT lateral distanoo from panel oontronl to oontor ol ^rav it \ . m 

Ai* low or surfaoo t* mums uppor surfaoo 0 

P p 1 1 p 

V offoottv o twisi audio, il#u 


Subscript* 


lU'l'Y contribution ot the l*od\ 

« dorivativ o with rospoot to angle ol attack, por dod 

I* dorivatn o with respect to angle of yaw . por deg 

b t derivative xv xtlx rospoot tv' aileron doflootton, por do*: 

derivative w »t lx rospoot to olovtitor deflection, por dog 

b ( . derivalix o xv it lx rospoot to rudder deflection, por deg 


Ml S' HIP I'U'N Ol blNO HOMY PKOOUXM 


Operation of tho bing llody program rexpiiroa that a wing »'r t*''»t\ configura 
lion, or both, bo sp«vtfi»'il tit torms ot trapoconlal constant prossuro panels I low 
conditions. pninnrilv audio of attack and Maoh number. aro also spooiln\t Ising 
thoso inputs, tho program roprosonts tin' wing as a surfaoo distribution of source*, 
and vorfioitx Mini tho Ixulv as a lino distribution of souroos and doublets ui a 
linearized potonltMl flow fiold . I'lto h'*vl\ is simulatO'l k'no'' m axial flow and onoo in 
fU'w produced h\ tho presence of tho wind l ho program compute* the poti'ntial. 
stoa-lx stalo prossuro oooffioiont on oaolt body panol and on tho uppor amt lower 
surfaoo of o. toll wind panol. t uid tin* dit'forontinl prossuro oooffioiont . \v . wluoh 

is tho difforonoo between fin* i v duos for tho U*wor ami uppor surfacca, tho total 

airoraft fotvo ami un'iin*nt coefficients aro obtaiii»*»t 

I'ln* bind body rosults aro based on tho potential flow equation. If subsonn' 
flow is desired, rosults aro obtainod bx tin* use of a coordinate -.\ stom transfer 
matton basod on tho I'ramltl illauort faotor. which produces an otpiivalont conftg 
urattoil m iiK'omprossiblo flow If auporsomo fU'w is required. a traiisform«*'l 
vviifiduration at a Maoh numbor of x i*. obteiuod. wluoh i*. cquixalcut tv' tho »l»'sir'St 
configuration and Maoh nunibor l'ho result? aro llion adjusted t»' tho desired Maoh 
numbor l so of tin* prodram in fho transoun* r«*dion is not valid because »*t \ i scons 
offivts that position tin* omboddod sln'ok wavos l'ho suggested range *'f applioation 
is thoroforo subsonio »*r supersonic t'l (I to 0 S:* or M l ;* tv* .1 0> bind body 
doos not aooount for uonlitioar flow oonditions suoli as \''rt»*\ formation or flow 


!.>«♦* 1 M l l l I * ! - 1—1 » Ljl : i i l— m 


separation and stall. riiiToforo. ivoffloli’iilii ilorlvrtl lining Wing IloU) will always 
be linear with respect to angle of ntUok. regardless of the type of flow actuallv 
experienced by the aircraft . 

The severest limitation of tins version of the Woodward t'armichsel method 
has been the number of panels that can be used. Wing surface description is 
limited to 100 panels Other version* accept the use of more panels but generally 
require substantially more computer time For symmetric configurations. Wing 
Hoily allows tin* entire unmbei of panels to be used to represent one half of the 
configuration . The resulting coefficients are then doubled . 

Hie leading and trailing edges of each panel may be swept, however, the 
two side edges must be parallel to the fret' stream Triangular panels are accept 
able if a atrcamwlse side is set equal to zero Dihedral . camber, and thickness 
may also be specified tor each panel. I'amber inputs are ctpilvaleut to the local 
angle of attack of each panel. When describing a computer model, it is important 
that all wing surfaces be subdivided into panels at the same span stations so that 
the panels along am streamwise row will have the same width and their stream 
wise side edges will coincide 

The program incorporates several features that reduce the Input preparation 
time. Data are arranged in namelist form for which most variables have default 
values Txeept for angular values, no specific unit sv stem is required because 
all results are uondiittcusioiialiccd \uother feature automatically di\ ides wing 
surfaces into the desinsl number of spanwtsc and chordwtse panels The paneling 
of the fuselage, including forebod) and afterlsHiv closures, is performed internally 
V plotting routine is available that allows verification of the paneling scheme input 
by the user 

On the i'lU' Cyber .*0, Model .M. Wing body iwithout overlay' requires 
1 77 ,000 octal words of central memory . ami most cases take less than three minutes 
of central processor time 


i OM U'd K V TU'\ S I Mill I 11 \\ M \ SIS 


Static stability derivativ es of flight \ chicles can be e\ aluated by Wing Ihniy 
with respect to angle of attack, sideslip, or control surface deflection* The 
following set of stability and control derivatives can be accurately determined 


Drag coefficients are 


1 l * l 'm ’ 1 l * 1 y‘ • v n • 1 n ’ 1 y * 1 v’ • * ml 1 y 
a a \s \s "ft u fl 0 \s \s 

e a a r r r r 

only first order approximations and are not suitable to evaluate vehicle performance. 

When appropriate, as in the case of oblique w ing vehicles, terms such as rolling 

moment due to elevator. C . , or pitching moment due to aileron, r . can 

v fS m tS 

e a 


be obtained 
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The F Ill transonic aircraft technology (TACT) research airplane waa analys'd 
by tlu> iib'o of (la* method. Stability changes reaulting from wing aw cop angle and 
Mach number variation* were determined anti the results compared favorably with 
the flight data included in reference 3. 

Five computer passes and three computer models were required to obtain a 
complete set of derivatives. To obtain derivatives with reaoect to angle of attack, 
sideslip, elevator deflection, aileron deflection, and rudder deflection , a separate 
computer pass was required for each case. The three models are shown in figure I 
and represent only the TACT aircraft’s planform. Wing thickness and camber do 
not affect the desired results. T'.e paneling scheme was purposely arranged such 
that groups of panels corrcapoi .1 in location and sue to control surfaces. Camber 
was then added to these control surface panels ami the resulting incremental changes 
in aircraft forces and moments indicate the control surface effectiveness . 

The program tends to overpredict the effect of camber inputs near the trailing 
edge. Therefore, when analyzing conventional trailing edge control surfaces, the 
control derivatives obtained are usually multiplied by a scale factor of 2/3. No 
scale factor is required when the control surface under consideration consists of an 
all movable wing surface. The TACT aircraft has both types of control surfaces, 
the rudder being a conventional trailing edge device, with the computed results 
being adjusted accordingly . and the aileron and elevator control being pro* ided by 
the ail moving horizontal stabilizer which required no adjustment. 

The first model (fig. 1(a)) was used to determine the longitudinal derivatives. 
The symmetric configuration option was employed, and, therefore, all the available 
panels could be used to represent one half of the aircraft . Vertical surfaces were 
not modeled . 

The second model (fig. 1(b)) was used te estimate lateral derivatives. The 
asymmetric optic n was used here si' that differential camber inputs could be made to 
the aileron panels, because of this, the paneling was much coarser for this model. 
Only derivatives with respect to aileron deflection arc obtained from this model. 

To obtain sideslip and rudder control derivatives, a third model was used. The 
vertical surface was described in the \ Y plane, as shown in figure 1(c). Angular 
changes in the X 7. plane, which would usually be considered angle of attack, repre 
sent sideslip to this model. The asymmetric option was used . 

Die horizontal surfaces, as well as the fuselage, were not included because 
their effect on the directional derivatives is small. On some aircraft the fuselage 
should be included in the determination of the directional derivatives, although on 
the FACT aircraft the fuselage is' distributed evenly about the center of gravity, 
because the coordinate axis system has been rotated 90° for ihis model, forces in 
the 7. direction labeled (* correspond to side forces on the vertical tail, anti momenta 

labeled (' m represent yawing moments. Using this analogy, sideslip and rudder 

control derivatives can be obtained. 

Some of the TACT stability and control derivatives determined with the Wing 
body program are shown in figure 2. Corresponding flight data are also shown. 


i i i 


The result* are for Mach l).‘t flow with the wing* awept to 2tf°. In general, the 
correlation between Winn hotly and flight i| u ta is good ■ 

The longitudinal derivativea In figure 1(a) demonstrate good agreement up to 
an angle of attack of 7°. Flow separation in expected above thin angle of attack. The 
predicted elevator control derivatives in figure 2(b) are generally higher than 
experimental value* . A contribution to thin derivative that is not accounted for by 
Wing hotly In the decrease in dynamic pressure in the wing wake. 

The aileron control derivatives in figure 2(c) are computed with respect to 
differential deflection of the rolling tail surfaces. Yawing moment due to aileron 
deflection, (' , from both flight and computed data is very small. The Wing 

"ft 

a 

Hotly results are computed neglecting the vertical tail. For this reason, flight data 
yield slightly proverse yaw while Wing Hotly predicts slightly adverse yaw. Results 
for rolling moment due to aileron deflection, C fl , correlate well. 

ft 

a 

The predicted yaw and roll moment derivatives with respect to rudder deflection 
(fig. 2(d)) al to show gootl correlation with flight data. 

A similar stability and control analysis has been made of the ADI oblique wing 
research vehicle at various asymmetric wing sweep angles. The results were used 
to verify the design and procure the airplane. 

FLOW FIHI.D DFTFKMIN ATION 


Wing Hoily has been used to estimate the potential flow direction in front of an 
aircraft in subsonic flow and to identify nonpotential flow in flight data. Kxnmples 
of each of these studies will be discussed. 


The upwash angle created forward of ihc TACT wing was calculated using Wing 
body. The computer model includes the TACT aircraft planform as well as a small 
fictitious wing panel located a few centimeters forward of the TACT wing's leading 
edge. Hy making tin* fictitious panel sufficiently small, its effect on the flow field 
around tin* TACT airplane can be neglected. The AC value of '.his panel will be 

calculated with the effect of the TACT upwash field. A second computer pass was 
required to determine the lift characteristics of the fictitious panel without the 
influence of the TACT aircraft . The differential pressure coefficient of this panel 
is shown as a function of angle of attack for both of these computer runs in figure 3. 
The increase in AC^ due to the presence of the TACT planform is a function of the 

upwash angle. For any given AC^ value, the difference in angle of attack between 


the two curves is equal to the upwash angle. This technique has also been used to 
suggest angle of attack vane corrections for the highly maneuverable aircraft 
technology (HiMAT) and A *0 aircraft. 


ti 


Another u*e of lit* W ut|{ Uodx program has bx»on t»» tdx'iitifx nonpotontial flx»w 
when u in oxporimonlul xlata. |i\ obtaining a go, si pxitx'iittal tluvry preaaure 

solution using W inf IUhIv • nonlinear f.ow ofteets will be absent Comiturobk 
Main Inkon in flight may ttteludo vxirtox . separatist, or ollior (low offeota When 
flight mnt potential pressure value* are x'omparxHl, nouliitoarilios m tin* flight Mata 
htvoino etx*ar . 

I'h«' fornialton x't tho loading xslge vortox on tho Ming glovo rogion of the I'Al'l' 
nnvrnft wm identified by oompartng tho flight meaatiriHt amt oaloulattxt ixitontial 
pressure values. I'tio Wing UoMx modi ' uaod tor tins stmt) is show n in liguro 4 
I'lio distribution ol patiols was hoax ilv xSMtoxMit rated near tho xx nig root t»» impr»»vo 
tho aeouraoy of tho prossnro value* computed in this rogion. I vportmontal pros 
anro Mata wore ohtanioM from tho full sonlo I'Al'l’ airoraft Prossnro onfioos were 
Ixvatisl m tho xx nig glox o rogion, amt inoasnromont s wor* taken primarily as a 
fnnotion ot .niglo of attack. 

Pho rosnlts from flight amt from tho px'tontial solution aro oompurxxt in fignro > 

Vt h*xx anglos of attaok it is obx ions that tho potontial thoory prox nloa a oloso 
apprx'ximation of tho monsuroM prossnro x alnos \s tho anglo of attaok is moroa soM 
a!x'\ o <P*. somo nonlmoar ofuvts ill tho flight xlata btvx'iito ovtxtent. I vpononoo 
suggosts tho formation of a x ortox at this ox'intltion . xvlnoh proMnoos tho nonlinear 
Mix orgonoo from tho potential solution 

t \l M U I'N N Wtll'S 

l sing tho Wing Uodx technique, an analytioal at mix of tho proiHxsxsl \ Mi' 
National Hypersonic Plight Research Facility lannoh dynamics was performed I'lio 
\ .' U‘ aircraft was intonMoM to Ih» oarrioM to a snitablo lannoh altitmlo hx tho 
|i :x.’ airplano amt tlion roloasvHt Pho mtorforonoo effxvt* of tho It b.’ airplano on 
tho \ Mi’ stability dorivattvoa Muring snoli a i'roposoil lannoh wore calculatxvl 
I’lningoi in pitolnng niomont . rolling momont . anM lift as a fnnotion of tho x ortioal 
soparatton t'otxvoon tho two vehicles was of primary intorost \lso. as a part of this 
study . changes in tho \ Ml' Mosign . pitch angle, amt launch pylon position xvero 
analx.ovt. i'hangos in control effectiveness were not MolorminiHt, lu*xxovor, those, 
as well as tho effects of \ Mi’ roll attitude. oonlM also have Ivon stmt tod using the 
Wing ItoMy program. 

Varti'Us \ ’ 4i' designs xvero modeled xxtth both wing ami boMx panels t’o 
insure that an adequate computci moxtel ItaM been obtained, tho \ Ml' airoraft xvas 
first annly . oM without the influence of the It carrier aircraft Wind tunnel data 
wore available for one design, and they aro compared with Wing UoMx results in 
figure tx. x i»s'«l correlation is shown bx'txvoon tho experimental ami thoorotioal values 
of longitudinal xMoftioionts l it. Wing Itxwty x*omputer moxtel of tho \ Mi’ airoraft was 
then usoxl as part of tho mated configuration model suailar to that ahoxvn in figure 7. 

Itxvauso this pri'gram allows only otto airoraft fiiaelage t»» bo uixxxlelexl . the 
11 airoraft was roprosoutxxl by a xving aurfaoo only \ separate sot x>f ox*mputer 
passes xvas made to xlotormino tho oamber for this xx tug surface so that tho model 
woulxl inoluxto the offxvt of tho H S7 fuselage on tin" lift distribution. First, tho 
It v.* airoraft was modeled xxtth both x\ nig and bodx panels I'lio reuniting Ax' ^ 
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values on the right wing were used as the input for a subsequent computer pass. 

A design option, which is available from Wing Body , was used to compute the 
camber values required to produce a given AC distribution. For this model. 

the fuselage and left whig were eliminated. Using the camber values generated 
by this computer pass, the right wing panels produced the same loading both 
spanwise and chordwise as the complete It 52 configuration. 

The 11-52 and X 24C models were then combined geometrically (fig. 7). Again 
it was necessary to cause the streamwise panel edges to aline at common span stations. 
The configuration in figure 7 shows the X 24C aircraft in the mated position with 
only 0.3 meter of clearance between its fuselage and the wing of (he It 52 aircraft. The 
X 24C fuselage centerline had to remain in the 7. 0 plane due to the program's input 

arrangement. Therefore, to simulate an increase in vertical separation distance 
between the two vehicles, the /. coordinates of the It 52 wing had to be increased. 

To obtain the effects of angle of attack on the X 24C aircraft, the angle of attaek 
of the B 52 aircraft had to remain constant. When the angle of attack input parameter 
is used, the angle of attack of both the X-24C and H 52 models is increased. To 
maintain a constant angle of attack for the B-52 wing, the camber values for the wing 
panels had to be reduced by the same amount as the angle of attack input parameter. 

To obtain the X 24C stability terms, it was necessary to sum the panel forces 
externally to the Wing Body program. The summation of panel forces performed 
automatically by Wing Body would include the lift and moments generated by the 
B 52 wing. Therefore, the following equations were used to total the X 24C aircraft's 
lift, roll, and pitch coefficients . 


C L = “V 


1 K 

♦ i V (AC ) A 
, h *-*, p n 

( n-1 K n 


BODY 


C ‘ (C «> V YC 'n 

BODY n-1 * n 


1 K 

C = (C ) ♦ — y (AC ) A (XC) 

m m BODY SS n=l *> n n n 


where K is the number of X 24C panels. A brief FOR TRAN computer program was 
used to perform this summation using punched cards of AC values generated 

directly by Wing Body . 

Results obtained using the configuration shown in figure 7 are given in figure 8. 
Data for two launch separation distances are shown: the mated position and the 
position after the X 24C aircraft has descended half the distance of the B 52 wing 
chord. The influence of the B 52 aircraft on the X-24C aircraft drops off quickly as 
the separating distance increases. At higher separation distances, the characteristics 
of t tie X 24C aircraft resemble the isolated X 24C results. 
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Vs well as prey iduig those ty p«*» of data, other aspect* of a proposed \ .'4t' 
launch were studied I'ln* span station at which the \ .’4i* aircraft would be attached 
was vnriiHl, hihI tin* romilli* are »tu«wn tn figure •» Only small difference* 

M CIV Observ »'»1 
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I'he force* and moments created bx trailing wake vortexes own be predicted using 
tin* \vnt|! body program. Velocity component* in tin* atmosphere that are normal to 
tin* fligtitpath of an aircraft own In* simulated tn tin* computer model of tin* \ etuelc 
I'ln* value ot a \ crtical velocity »*oinpoiu*nt whon related to the airspeed o r the trailing 
\ chicle determines an additional increment of angle * * attack at the kVrrriiHiiHlmn 
region of the wind Vito \ «*t*tn*al velocity . which varies across the span of the 
trailing aircraft, can In* effectively approx minted t»\ modeling wing twist tn tin 
trailing aircraft I'ln* value of tin* induced twist angle. v . ia deny ed as follows 

y tan ' (\ Vi 

V wake vortex induct'd angle of attack can he calculated to** each s»rcamwise row ot 
panels ot a trailing aircraft computer model rhereforc. the computed trailing 
aircraft force*. amt momenta will inchnte the effects ot the \ ertical airflow 

rite rolling nnmient . . er *ated on a trailing aitvraft m the wake of a M 74* jet 

transport was predicted using this method I'ln* effect i\ e twist of the trailing air 
craft computer model was determined from a velocitv profile of the U 7 47 wake 
vortex***, which were obtained from flight experiments Vhla H 4 7 wake vortex 
velocity profile is slnnvtt in figure 10ya> 

The rolling moment created on a trailing aircraft fixing tn such a flow condition 
varies with Us location tn the vortex, however, the maximum rolling moment is 
created when the trailing aircraft’s centerline is positioned m the middle ot the vortex 
In addition, the rolling moment produced on a trailing aircraft is a function of the 
platiform of the trailing aircraft. 

I'ln* effects of tin* It .'47 wake vortex flow field on ls>th a I’ .17 and I cartel air 
craft were analy zed tn this study . i* output er models **1 the two aircraft consisted ot 
only w ing panels I'ln* fuselages were neglected because of their smr.il Contribution 
to the rolling moment coefficients Hie horizontal tail surface, were tv.t tinMclcd 
because the flow field around these surfaces is dominat'd by the townwaah of the 
wing rite T .17 wing is generally rectangular, while the 1 earjet wing has a 
t* 7> percent greater span, a 17*. V’ leading edge sweep, and a taper ratio of 0 :*l 
Moth computer models were divided spanwise Into It) row s of panels For each 
streanwise row of panels. • vortex induced angle of attack was calculated t*a>«xl on 
the It 7 4 7 flight data shown tn figure lOyai and a trailing aircraft velocity . V . of 
*h! in sec The eff.vtix e twist angles produced t*y the peak vortex \ eUvities were 
a Unit 10° 


The results are given In figure 10(b), The rolling moment of the T 37 nirernft 
dropped off quickly as the aircraft moved away from the vortex core. AIho, the 
maximum rolling moment created on the Learjet aircraft waa lower than that created 
on the T 37 aircraft . 

Thia atudy allowed an analytical comparison of two aircraft planforma subjected 
to an identical vortex flow field. In flight, it would be impoaaible to recreate a 
specific velocity profile for different trailing aircraft. Alao. it would not be aafe to 
subject some aircraft to such a wake vortex upset experimentally. 


CONCLUDING REMARKS 


Through experience gained at DKRC in the application of the Wing body aero 
dynamic analysis program, several investigations of aircraft characteristics have 
been made analytically Examples of the following studies have been presented: 

1 . Determination of aircraft stability derivatives with respect to angle of attack, 
angle of sideslip, and control surface deflection . 

2. Determination of the upwash angle in front of an aircraft. 

3. Identification of the effects of nonpotentiai flow in experimental data. 

4. Determination of interference effects of a carrier aircraft on an air 
launched vehicle. 

5. Determination of the loading produced on an aircraft by wake vortex flow. 

Rased on comparisons with experimental data presented in this report. Wing 
Di ly analytical results may be considered reliable. Results have been used to plan 
wind tunnel or flight testing, or as verification of experimental work. Some 
situations that cannot be investigated readily or safely by experiment can be studied 
analytically using W’ng body. 


Drytien Flight Research Center 

National Aeronautics and Space Aiimmsfration 
Edwards, Calif., October -f I 1 ?# 
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Figure 2. Concluded . 
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Figure 3. Pressure coefficient versus angle of attack of fictitious panel. 


18 



-.14 


Wing -Body 
Flight 


0^ measured 

at an upper 
surface 
pressure 
orifice 
near the 
leading edge 






0 2 4 6 8 10 12 

Angle of attack, deg 


j 


Figur e 5. Comparison of potential theory results and flight data. 
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higure 8. ( omparison of X 2-k' characterisin' s in pretence of H 52 aircraft. 
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Fiyure S. Concluded. 
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